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METHOD AND SYSTEM FOR INDICATING ANISOTROPIC RESISTIVITY IN 

AN EARTH FORMATION 

1. BACKGROUND OF THE INVENTION 

1.1 Field of the Invention 

The present invention is generally directed to the analysis of underground earth 
formations, and, more particularly, to the determination of formation resistivity properties 
and/or profiles. 

12 Description of Related Art 

Knowledge of the resistivity properties of underground earth formations is 
fundamental in the evaluation and characterization of potential and existing hydrocarbon- 
bearing reservoirs. The terms "logging" and "log" are respectively used herein to 
designate any sequential measurements and the recording of such measurements as a 
function of depth and/or time of variations in a given characteristic of the formations 
around a borehole. 

Formation resistivity anisotropy is a characteristic of reservoirs that can complicate 
formation evaluation. Many reservoir rocks exhibit resistivity anisotropy, especially when 
saturated with oil. There are several mechanisms, which can produce this anisotropy, 
among which are very thin sand-shale laminations, depositional changes in clean 
sandstone, and wind-distributed sands (aeolian formations). See Rubin, D. M., Cross 
bedding, bedforms t and paleocwrrents, SOCIETY OF ECONOMIC PALEONTOLOGISTS AND 
Mineralogists, Concepts in Sedimentology and Paleontology, 1; Klein et al., The 
petrophysics of electrically anisotropic reservoirs, Transactions of the SPWLA TfflRTY- 
Sdcth Annual Logging Symposium, Paris, France, June 26-29, 1995, paper HH. 
Resistivity measurements from within a borehole drilled through these formations have 
been made by a number of techniques employing the use of well-known instruments, 
including induction and laterolog (electrode) type devices. 

In induction logging, a transmitter coil mounted on a well tool is lowered into a 
borehole and energized by an alternating current. Indications are then obtained of the 



BNSOOCID: <QB 23fl891iA_l_> 



influence of sounding formations on the electromagnetic Held established by the coil 
UsuaJly such indications are obtained by observing the voltage induced in a receiver coil 
also mounted on the well tool in coaxial relationship with the transmitter coil and axially 
spaced apart therefrom. Conventional induction well tools include at least one transmitter 
and a receiver coil array. Such an induction tool is illustrated and described in U.S. Pat. 
No. 5,157,605, assigned to the present assignee. 

The laterolog technique comprises making a measurement with a well tool having a 
current emitting electrode and voltage electrodes one either side thereof which are 
operated to force current into the formation. The current is focused into the formations 
and returns locally, i.e., to the tool body. The formation resistivity is determined by 
measuring the current flow. Conventional laterolog tools are illustrated and described in 
U.S. Pats. Nos. 5,754,050, 5,852,363, 6,023,168, and Documents EP-478409, EP-544583 
all assigned to the present assignee. Document EP-544583 describes a laterolog tool 
including both a "deep" investigation mode and a "shallow" (LLS) investigation mode. 

The sensitivity of a resistivity logging tool to resistivity anisotropy depends on the 
Physics of the measurement and the relative dip angle of the well bore of the formation 
For example, 2MHz induction-type tools have no sensitivity to resistivity anisotropy when 
the formation layers are perpendicular to the well bore, but will be quite sensitive at 
relative dip angles above 60 degrees. See Luting et al., Processing and modeling 2 MHz 
resistivity tools in dipping, laminated, anisotropic formations, Transactions of the SPWLA 
Thirty-Fifth annual Logging Symposium, Tulsa, Oklahoma, June 19-22, 1994 paper 
QQ. 

Induction logging tools are insensitive to resistivity anisotropy at low relative dip. 
The current density set up in the formation is limited to the horizontal direction. See 
Moran et al., Effects of formation anisotropy on resistivity-logging measurements, 
Geophysics, V. 44, No. 7, July 1979, 1266-86; Anderson et al., The Response of Induction 
Tools to Dipping Anisotropic Formations, Transactions of the SPWLA Thirty-Sdcih 
Annual Logging Symposium, Paris, France, June 26-29, 1995, paper D. At higher 
Native dips the measured resistivity will combine both vertical and horizontal 
components. However, the effect does not depend on the coil array spacing, so the 
components are not resolvable. See Anderson et al., The Response of Induction Tools to 
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Dipping Anisotropic Formations, Transactions of the SPWLA Thirty-Sdcth Annual 
Logging Symposium, Paris, France, June 26-29. 

Laterolog logging tools have some sensitivity to the vertical component of 
resistivity. See Moran et al., Effects of formation anisotropy on resistivity-logging 
measurements, GEOPHYSICS, V. 44, No. 7, July 1979, 1266-86; Chemali et al., The effect 
of shale anisotropy on focused resistivity devices, Transactions of the SPWLA TWENTY- 
Eighth Annual Logging Symposium, London, England, June 30-July 2, 1987, paper H. 
However, with the dual-laterolog, the deep laterolog was not very sensitive to anisotropy. 
The shallow laterolog, which returns current to nearby electrodes, is more sensitive to 
anisotropy (See Moran et al., 1979), but invasion will reduce the effect of anisotropy. 

It is desirable to obtain a reliable method and system for determining the resistivity 
properties and profiles of potential hydrocarbon-bearing zones in subsurface formations. 
Thus, there remains a need for improved methods and techniques to indicate formation 
resistivity anisotropy and evaluate the vertical component of resistivity. 

2. SUMMARY OF THE INVENTION 

Methods and systems are disclosed for combining induction tool and laterolog tool 
measurements to indicate that the formation resistivity is anisotropic and estimate the 
vertical component of resistivity even in the presence of invasion. 

One aspect of the invention provides a method for determining a characteristic of a 
subsurface earth formation surrounding a borehole. The method comprises calculating at 
least one parameter representative of a property of the formation using log data obtained 
with a first well tool disposed within the borehole; modeling log data from the at least one 
calculated parameter; and comparing log data obtained with a second well tool disposed 
within the borehole against the modeled log data to determine the formation characteristic. 

Another aspect of the formation provides a method for determining a characteristic 
of a subsurface earth formation surrounding a borehole. The method comprises disposing 
first and second well tools within the borehole to obtain measurement dat a; producing at 
least one parameter representative of a property of the formation from the measurement 
data obtained with the first well tool; producing theoretical data from the at least one 
produced parameter; and comparing the theoretical data against the measurement data 

obtained with the second well tool to determine the formation characteristic. 

•i 
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Another aspect of the invention provides a well logging system including first and 
second well tools adapted to be moveable through a borehole and an apparatus adapted to 
be coupled to the well tools, the apparatus being adapted to respond to data obtained with 
the well tools to determine a characteristic of a formation. The apparatus comprises means 
for calculating at least one parameter representative of a property of the formation using 
log data obtained with the first weU tool disposed within the borehole; means for modeling 
log data from the at least one calculated parameter; and means for comparing log data 
obtained with the second well tool disposed within the borehole against the modeled log 
data to determine the formation characteristic. 

Yet another aspect of the invention provides a program storage device readable by 
a processor and encoding a program of instructions including instructions for performing 
operations in accord with the invention. 

3. BRIEF DESCRIPTION OF THE DRAWINGS 

Other aspects and advantages of the invention will become apparent upon reading 
the following detailed description and upon reference to the drawings in which: 

Figure 1 is a schematic diagram of the lateroiog energization patterns obtained with 
a conventional lateroiog tool. 

Figure 2 is a schematic diagram of the current lines corresponding to the lateroiog 
patterns of Figure 1. 

Figure 3 is a graph showing the response of the five lateroiog arrays of Figures 1 
and 2 in a homogeneous anisotropic formation. 

Figures 4-7 are graphical illustrations of the responses of the five lateroiog arrays 
of Figures 1 and 2 in invaded, anisotropic formations. 

Figure 8 is a graph showing the comparison of a lateroiog tool response in an 
anisotropic formation vs. its response in an isotropic formation in accord with the 
invention. 

Figure 9 illustrates several plots of a ratio used in the calculations of the invention. 
Figure 10 illustrates the ratio plots of Figure 9 with an applied scaling factor. 
Figure 11 shows a crossplot for seventy ratios pertaining to cases of formation 
parameters in accord with the invention. 
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Figure 12 shows an average value of a parameter for the ratios of Figure 11 in 
accord with the invention. 

Figure 13 is a graph showing a computed formation vertical resistivity component 
in accord with the invention. 

Figure 14 is a graph showing computed induction tool and laterolog tool logs in 
accord with the invention. 

Figure 15 is a graph showing the inversion of the induction tool data of Figure 14 
using a laterolog model in accord with the invention. 

Figure 16 is a graph showing the vertical resistivity component derived from the 
data of Figures 14 and 15 in accord with the invention. 

Figure 17 shows field logs obtained from an induction tool and a laterolog tool in 
accord with the invention. 

Figure 18 shows the vertical resistivity component derived from the data of Figure 
17 in accord with the invention. 

Figure 19 is a schematic diagram of a well logging system in accord with the 
invention. 

Figure 20 is a flow diagram of a method in accord with the invention. 
Figure 21 is a flow diagram of another method in accord with the invention. 

4. DETAILED DESCRIPTION OF SPECIFIC EMBODIMENTS 

In the interest of clarity, not all features of actual implementation are described in 
this specification. It will be appreciated that although the development of any such actual 
implementation might be complex and time-consuming, it would nevertheless be a routine 
undertaking for those of ordinary skill in the art having the benefit of this disclosure. 

Figure 1 shows the energization patterns, referred to as "laterologs", obtainable 
from a laterolog-type tool 10 such as those disclosed in U.S. Pats. Nos. 5,754,050, 
5,852,363, 6,023,168, and Document EP-478409. The deepest array, mode 5, is 
considerably deeper than the LLS and should have considerable sensitivity to anisotropy. 
Figure 2 shows the current paths on several of the arrays of the tool 10. Note that the 
current flow lines turn parallel to the tool 10 axis and return locally, thereby increasing the 
sensitivity to the vertical resistivity component. 



Figure 3 shows the response of the tool 10 arrays (1-5) to an anisotropic 
homogeneous formation. The ratio of the apparent resistivity Rapp over the formation 
horizontal resistivity component Rt„ is plotted against the ratio of the formation vertical 
resistivity component Rt v over Rt„. This separation shows that the shallow arrays are 
much more sensitive to an anisotropic formation than the deeper arrays. Although the 
homogeneous medium anisotropy response looks favorable, what really matters is the 
response in the presence of invasion. 

Figures 4 through 7 show four representative cases of formation resistivity 
measurements obtained with a laterolog tool in an anisotropic formation having invasion 
out to a radius of 60 inches. These figures show ratios of the formation horizontal 
resistivity component Rth over the formation invasion resistivity Rxo of 0.1, 0.2, 5, and 
10. For all the cases Rtv = 10 x Rth. The invaded zone is isotropic with resistivity Rxo. 
In Figures 4 and 5, Rxo > Rth. In Figures 6 and 7, Rxo < Rth. For all these cases, the 
invasion response confuses the anisotropy response. 

Figure 8 shows the anisotropic case of Figure 7 superimposed with the laterolog 
tool response to an isotropic virgin zone with formation resistivity Rt equal to Rth. As 
seen in Figure 8, the two curves representing Array 5 (RLA5) are almost parallel over the 
range of invasion radius out to 60 inches. The other shallower isotropic-vs.-anisotropic 
cases show separation at shallow invasion, but little or no separation at deeper invasion 
radii. The reasonable constancy of the difference between the isotropic and anisotropic 
Array 5 response suggests that the ratio of anisotropic to isotropic response or other ratios 
provide useful information. 

Designating Raa as the laterolog Array 5 log in an anisotropic formation, and Rai 
as the laterolog Array 5 log reading in an isotropic formation, the following ratio is 



derived: 

Rai 



(Raa -Rai) 



Figure 9 shows a plot of the above ratio for Rxo =1 and Rth =10 for A values from 1.4 to 
10, where 



X = 




Rtv 
Rth 
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It should be noted that the spread is almost proportional to X . This suggests that 
the ratio 

XRai 
Raa - Rai 

should be nearly constant. This ratio is shown in Figure 10. 



If the ratio 

A= ^ (1) 
(Raa -Rai) 

is nearly constant for this case, then other values of Rxo and Rth should be looked at. 

Figure 11 shows a crossplot of 7 combinations of Rxo and Rth. There are 10 values of A 

ranging from 1.4 to 10. The spread in A is remarkably small. Figure 12 shows the average 

value of A from the data of Figure 1 1 . If an average value of A is taken to be a constant (or 

fitted function of the radius of invasion ri), then the unknown X can be solved for 

XsiA (Raa-Rai) (2) 
Rai 

One unknown remains in equation (2) - Rai. 

As mentioned above, induction tools have litde or no sensitivity to Rtv in a vertical 
well. In such a case, the induction tool log data is inverted for Rxo, Rth, and ri using a 
lookup table, computer algorithm, model, or any other inversion means as known in the 
art. Also known in the art is the use of modeling techniques to simulate specific tool 
responses from given data. By applying a one-dimensional (ID) model of the laterolog 
tool response to the formation parameters obtained from the induction tool, Rxo, Rth, and 
ri, the laterolog log data in an isotropic formation (Rai) can be inferred. 

Then equation (2) can be solved for X , and then 

Rtv=A 2 Rth. (3) 

Figure 13 shows such a solution using a fitted value of A as a function of ri. The values 

for A may be precomputed. The Value for A ID model for the laterolog tool was derived 
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using a simple average ./-function, which is only a gross estimate for an example The 
estimate of Rtv is within a factor of 2 of the actual value of Rtv in the uninvaded zone 
Although a factor of 2 may be considered inaccurate for estimates of petrophysical 
quantities, such accuracy is ample for determining that anisotropy exists in a vertical well 
The accuracy may be improved by refining A. Both induction and laterolog tools have 
effects in 2D formations (layers and invasion in the layers) that are not accurately modeled 
with ID forward models. To test the usefulness of the ID models used in the recreation of 
Rai. it is advisable to try the algorithm in a 2D-formation model. 

Figure 14 shows computed induction and laterolog tool logs in a simple 2D 
formation with 3 invaded, anisotropic layers of varying thickness. These layers have 
identical values of Rtv, Rth, Rxo, and ri. The thickness varies from 50 ft. to 10 ft The 
induction tool log data reflects the receiver to transmitter spacings (10 in, 90 in.) for the 
particular tool. 

Figure 15 shows a modeled laterolog tool log derived from the induction tool log 
data of Figure 14. The modeled log data represents the theoretical log values obtainable 
with the laterolog tool disposed within the isotropic formation. The modeled laterolog 
log data are then compared against the actual or apparent laterolog log data to compute 
Rtv as described above. Figure 16 shows the resulting Rtv log computed from the data of 
Figures 14 and 15. 

The Rtv log of Figure 16 shows anomalies in the non-invaded isotropic beds 
between the invaded anisotropic beds. As seen in Figure 14, the laterolog tool logs 
display shoulder effect in these beds and on either side of the first and third beds Any 
source of error that leads to an increased separation between the laterelog tool logs 
predicted from the induction tool parameters and apparent laterolog tool logs will be 
interpreted as anisotropy. Overall, however, intelligent interpretation should remove 
most ambiguity. 

A field test was conducted comparing a laterolog tool against an induction tool in a 
formation with expected resistivity anisotropy. Figure 17 shows the induction tool and 
laterolog tool logs in a pay sand. Figure 18 shows the indicated vertical resistivity Rtv 
along with the induction tool Rth, Rxo, and ri inversion in a zone that exhibits the 
anisotropy signature. It also exhibits shoulder effect just as in the computed example of 
Figure 16. 
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4.1 Program Storage Device 

It will be apparent to those of ordinary skill having the benefit of this disclosure 
that the present invention may be implemented by programming one or more suitable 
general-purpose computers having appropriate hardware. The programming may be 
accomplished through the use of one or more program storage devices readable by the 
computer processor and encoding one or more programs of instructions executable by the 
computer for performing the operations described above. The program storage device may 
take the form of, e.g., one or more floppy disks; a CD ROM or other optical disk; a 
magnetic tape; a read-only memory chip (ROM); and other forms of the kind well-known 
in the art or subsequendy developed. The program of instructions may be "object code," 
i.e., in binary form that is executable more-or-less directly by the computer, in "source 
code" that requires compilation or interpretation before execution; or in some intermediate 
form such as partially compiled code. The precise forms of the program storage device 
and of the encoding of instructions are immaterial here. 

4.2 Logging System 

The present invention can be implemented in a logging system 200 with a first well 
tool 210 and a second well tool 215 adapted to be moveable through a borehole and an 
apparatus 220 coupled to the tools as shown in Figure 19. Certain conventional details are 
omitted in Figure 19 for clarity of illustration. The apparatus 220 comprises a surface 
computer (such as the general-purpose computer and program storage device described 
above) coupled to the well tools 210, 215 by a wireline cable 225 or linked to the tool 
string's surface instrumentation as known in the art. 

The first well tool 210 may be any laterolog-type tool that returns current locally, 
i.e., to the sonde body. The second well tool 215 may be any propagation or induction tool 
that provides multiple curves, preferably providing at least three independent 
measurements to solve for Rxo, Rt, and ri. In addition to the U.S. Patents mentioned 
above, other useable tools include the instruments described and illustrated in U.S. Pats. 
Nos. 5,339,037 and 4,968,940, both assigned to the present assignee. 

The formation parameters can be determined in real-time by sending the 
measurement data to the surface as they are acquired, or it can be determined from a 



recorded-mode by recording the data on a suitable recordable medium. As known in the 
art. the measured data are transmitted from the tools 210, 21.5 to the surface computer 220 
by electronics (not shown) housed in the tools. The data may be sent to the surface 
computer 220 along the wireline cable 225 or by mud-pulse telemetry if implemented in a 
loggmg-while-drilling configuration. I, will be understood that alternative means can be 
employed for communicating the acquired data to the surface as the precise form of 
communication is immaterial to the implementation of the disclosed techniques. 

Once received by the surface computer 220, the data can be recorded, processed or 
computed as desired by the user to generate the formation resistivity profile. The profile 
can then be recorded on a suitable output record medium. Alternatively, some or all of the 
processing can be performed downhole and the data can be recorded uphole, downhole or 
both. It will be understood that the present invention can be implemented in conjunction 
with any suitable technique for keeping track of the tool 210, 215 depth within a borehole. 

Figure 20 illustrates a flow diagram of a method for determining a characteristic of 
a subsurface formation surrounding a borehole in accord with the invention. The method 
comprises calculating at least one parameter representative of a property of the formation 
using log data obtained with a first well tool disposed within the borehole 300; modeling 
*> log data from the at least one calculated parameter 305; and comparing log data obtained 
with a second well tool disposed within the borehole against the modeled log data to 
determine the formation characteristic 3 10. 

Figure 21 illustrates a flow diagram of a method for determining a characteristic of 
a subsurface formation surrounding a borehole in accord with the invention. The method 
comprises disposing first and second weU tools within the borehole to obtain measurement 
data 400; producing at least one parameter representative of a property of the formation 
from the measurement data obtained with the first well tool 405; pnxhicing theoretical data 
from the at least one produced parameter 410; and comparing the produced theoretical data 
against the measurement data obtained with the second well tool to determine the 
formation characteristic 415. 

While the methods and apparatus of this invention have been described as specific 
embodiments, it will be apparent to those skilled in the art that variations may be applied 
to the structures and in the steps or in the sequence of steps of the methods described 
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herein without departing from the concept and scope of the invention. All such similar 
variations apparent to those skilled in the art are deemed to be within this 
concept and scope of the invention as defined by the appended claims. 
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CLAIMS 

A method for determining a characteristic of a subsurface earth formation 
surrounding a borehole, comprising: 

(a) calculating at least one parameter representative of a property of the 
formation using log data obtained with a first well tool disposed within 
the borehole; 

(b) modeling log data from the at least one calculated parameter; and 

(c) comparing log data obtained with a second well tool disposed within 
the borehole against the modeled log data to determine the formation 
characteristic. 

The method of claim 1. wherein step (b) comprises modeling log data 
representative of log data theoretically obtainable with the second well tool 
disposed within the borehole. 

The method of claim 1, wherein step (c) comprises calculating a ratio between 
the log data obtained with the second well tool and the modeled log data. 

The method of claim 1, wherein step (c) comprises calculating the difference 
between the log data obtained with the second well tool and the modeled log 
data. 

The method of claim 4, wherein step (c) comprises multiplying the calculated 
difference by a predetermined factor. 

The method of claim 1. wherein the first well tool is an induction-type tool and 
the second well tool is a laterolog-type tool. 

The method of claim 1, wherein the determined formation characteristic is a 
resistivity value or resistivity profile and the method is performed during or 
after drilling of the borehole. 
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8. A method for determining a characteristic of a subsurface earth formation 
surrounding a borehole, comprising: 

(a) disposing first and second well tools within the borehole to obtain 
measurement data; 

(b) producing at least one parameter representative of a property of the 
formation from the measurement data obtained with the first well tool; 

(c) producing theoretical data from the at least one produced parameter of 
step (b); and 

(d) comparing the theoretical data of step (c) against the measurement data 
obtained with the second well tool to determine the formation 
characteristic. 

9. The method of claim 8, wherein step (c) comprises producing theoretical data 
representative of measurement data obtainable with the second well tool 
disposed within the borehole. 

10. The method of claim 8, wherein step (d) comprises calculating a ratio between 
the measurement data obtained with the second well tool and the produced 
theoretical data. 

11. The method of claim 8, wherein step (d) comprises calculating the difference 
between the measurement data obtained with the second well tool and the 
produced theoretical data. 

12. The method of claim 1 1, wherein step (d) comprises multiplying the calculated 
difference by a predetermined factor. 

13. The method of claim 8, wherein the first well tool is an induction-type tool and 
the second well tool is a laterolog-type tool. 

14. The method of claim 8, wherein the determined formation characteristic is a 
resistivity value or resistivity profile and the method is performed during or 
after drilling of the borehole. 
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A well logging system including first and second well tools adapted to be 
moveable through a borehole and an apparatus adapted to be coupled to the 
well tools, the apparatus being adapted to respond to data obtained with the 
well tools to determine a characteristic of a formation, the apparatus 
comprising: 

means for calculating at least one parameter representative of a property of 
the formation using log data obtained with the first well tool disposed 
within the borehole; 

means for modeling log data from the at least one calculated parameter; 
and 

means for comparing log data obtained with the second well tool disposed 
within the borehole against the modeled log data to determine the 
formation characteristic. 

The system of claim 15, the modeling means further comprising means for 
modeling log data that is representative of log data theoretically obtainable 
with the second well tool disposed within the borehole. 

The system of claim 15, the comparison means further comprising means for 
calculating a ratio between the log data obtained with the second well tool and 
the modeled log data. 

The system of claim 15, the comparison means further comprising means for 
calculating the difference between the log data obtained with the second well 
tool and the modeled log data. 

The system of claim 15, wherein the first well tool is an induction-type tool 
and the second well tool is a laterolog-type tool. 

The system of claim 15, wherein the determined formation characteristic is a 
resistivity value or resistivity profile. 
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